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ABSTRACT 

Approximately 20-30% of 1 .4 ^ z ~ 2.5 galaxies with Xyega < 22 detected with Spitzer MIPS at 24 /im show 
excess mid-IR emission relative to that expected based on the rates of star formation measured from other 
multiwavelength data. These galaxies also display some near-IR excess in Spitzer IRAC data, with a spectral 
energy distribution peaking longward of 1.6/im in the rest frame, indicating the presence of warm-dust emis- 
sion usually absent in star forming galaxies. Stacking Chandra data for the mid-IR excess galaxies yields a 
significant hard X-ray detection at rest-frame energies > 6.2 keV. The stacked X-ray spectrum rises steeply 
at > 10 keV, suggesting that these sources host Compton-thick Active Galactic Nuclei (AGNs) with column 
densities A/h ~ 10 24 cm~ 2 and an average, unobscured X-ray luminosity L2_8keV ~ (1-4) x 10 43 erg s _1 . Their 
sky density (~ 3200 deg~ 2 ) and space density (~ 2.6 x 10~ 4 Mpc~ 3 ) are twice those of X-ray detected AGNs 
at z ~ 2, and much larger than those of previously-known Compton thick sources at similar redshifts. The 
mid-IR excess galaxies are part of the long sought-after population of distant heavily obscured AGNs predicted 
by synthesis models of the X-ray background. The fraction of mid-IR excess objects increases with galaxy 
mass, reaching ~ 50-60% forM ~ 10 M©, an effect likely connected with downsizing in galaxy formation. 
The ratio of the inferred black hole growth rate from these Compton-thick sources to the global star formation 
rate at z = 2 is similar to the mass ratio of black holes to stars in local spheroids, implying concurrent growth 
of both within the precursors of today's massive galaxies. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: active — X-rays: galaxies 



> 

m 
oo 

(N 

in 
o 
i> 
o 



X 



1. INTRODUCTION 

Active Galactic Nuclei (AGN) and galaxies are intimately 
connected. It is generally thought that virtually all nearby 
early type galaxies (and perhaps all relatively massive galax- 
ies) contain a supermassive central black hole (BH), indicat- 
ing that all massive galaxies may have experienced an AGN 
phase in their past (Kormendy & Richstone 1995; Magorrian 
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et al. 1998). Moreover, remarkably tight correlations have 
now been unveiled between central supermassive BH and the 
mass or velocity dispersion of the bulge in their host galaxies 
(Ferrarese & Merritt 2000; Gebhardt et al. 2000), with typi- 
cal mass ratio M B H/M bulge ~ (1-2) x 10~ 3 (e.g., McLure & 
Dunlop 2001; Ferrarese et al 2006). This suggests but does 
not establish that AGN and star formation activity may have 
been concurrent. The origin of these correlations is in fact not 
entirely clear, and has been the subject of active debate during 
recent years. It appears to be connected to the positive feed- 
back that the BH can exercise on its host galaxy, becoming ef- 
fective if the BH exceeds some critical mass limit and thus lu- 
minosity (Silk & Rees 1998; di Matteo, Springel & Hernquist 
2005; Springel, di Matteo & Hernquist 2005; Schawinski et 
al. 2006). If the energetic emission from accretion around the 
BH is large enough, the galaxy can be cleared out of its gas 
either via heating (e.g., Ciotti & Ostriker 1997; 2007) or ac- 
cretion related outflows as in BAL quasars (e.g., Chartas et al. 
2007), thereby preventing further gas accretion. 

In order to probe BH growth in the distant universe efficient 
tracers of AGN activity are needed. Obtaining a complete 
census of active galaxies is complicated due to AGN obscura- 
tion by circumnuclear material. During the formation epochs 
of galaxies, it can be well expected that a greater abundance 
of gas and its chaotic motions due to merging and interactions 
cause the typical AGN obscuration to be much larger than lo- 
cally. In that case, a substantial fraction of the cosmic BH 
growth history could be hidden from view and not yet taken 
into account. It is also currently unknown if the "Unification 
Paradigm" describing the AGN phenomenon (e.g., Antonucci 
1993; Urry & Padovani 1995) holds at high redshift, in par- 
ticular if the ratio of obscured to unobscured AGNs (or, more 
accurately, the distribution of obscuring column densities) is 
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the same as observed locally. More generally, the shape of the 
cosmic X-ray background, with its peak toward 30 keV (Mar- 
shall et al. 1980), suggests that much of the BH-accretion lu- 
minosity at substantial redshifts is obscured (e.g., Fabian & 
Iwasawa 1999). 

The deepest X-ray surveys available to date, performed 
with the Chandra satellite in the Chandra Deep Field 
South (Giacconi et al. 2002; GOODS-South field, GOODS- 
S hereafter), Chandra Deep Field North (Alexander et 
al. 2003a; GOODS-North, GOODS-N hereafter), and with 
XMM-Newton (e.g., Hasinger 2004 in the Lockman hole) re- 
solved most of the X-ray background at energies 0.5-6 keV, 
but a substantial fraction of the higher energy background re- 
mains unresolved (Krivonos et al. 2005; Worsley et al. 2005; 
2006; Hickox & Markevitch 2007). Synthesis models have 
been built (Comastri et al. 1995; Gilli et al. 2001 ; 2007) show- 
ing that in order to account for the high energy X-ray back- 
ground, where most of the accretion energy in the universe 
is being emitted, one has to postulate the existence of sub- 
stantial populations of high redshift heavily obscured Comp- 
ton thick AGNs (with column densities of A^h ^ 10 24 cm ). 
Such a population of Compton thick sources at high redshift 
still awaits discovery (e.g., Barger et al. 2007). Hard X-ray se- 
lection, the best means currently to find such sources at high 
redshift, is limited by telescope sensitivities and would allow 
one to find only the brightest sources and with non-extreme 
column densities (e.g., Brandt & Hasinger 2005). As an ex- 
ample, the archetype Compton thick AGN NGC 1068, with 
Lx > 10 44 erg s , would have a hard X-ray (2-8 keV) flux 
of only ~7x 10~ 18 erg s _1 cm~ 2 at z = 2, more than an or- 
der of magnitude fainter than that of sources detectable in the 
deepest (2 Ms) Chandra survey, in the GOODS-North field. 

On the other hand, some of the energy absorbed by the ob- 
scuring material is reprocessed and re-emitted at longer wave- 
lengths, providing the opportunity to detect these sources via 
their IR emission. For example, the nucleus of NGC 1068 
would have a flux density of 25/iJy at 24 /im if placed at z = 2 
— very faint but still detectable in the ultradeep Spitzer MIPS 
imaging taken as part of GOODS. One of the primary goals 
of GOODS was indeed the discovery of previously unknown 
populations of distant AGNs (e.g., Treister et al. 2004) by 
combining the deepest multiwavelength imaging brought to 
bear upon these fields. Several recent studies have focused 
on finding distant obscured AGNs by means of mid-IR or far- 
IR selection, in GOODS and other deep Spitzer surveys. For 
example, selection by means of red IRAC (3.5-8/im) colors 
were presented by Lacy et al. (2004), Stern et al. (2005), and 
similarly the selection by means of power-law spectral en- 
ergy distributions (SEDs) over IRAC bands (Alonso Herrero 
et al. 2006; Donley et al. 2007) have provided AGN candidate 
samples with only partial overlap with X-ray selected sources. 

The peak of QSO activity has been established to take place 
at z ~ 2 (e.g., Schmidt 1968; Croom et al. 2004; Hasinger 
et al. 2005). Therefore, it is quite natural to look for the 
growth of BH in massive galaxies at the same epoch. In- 
deed, the massive z ~ 2 galaxies selected at submillimeter 
wavelengths (SMGs) are undergoing star-formation activity 
and at the same time host luminous AGNs, just the signatures 
expected for joint stellar-BH growth (e.g., Alexander et al. 
2005a,b; Borys et al. 2005). However, these sources are rel- 
atively rare and appear to represent the most strongly star- 
forming galaxies at these epochs (e.g., Chapman et al. 2005) 
and it is important to see if this joint growth is also occurring 
in more typical massive galaxies at z ~ 2. 



This paper describes a source population serendipitously 
discovered while studying the multiwavelength emission 
properties of distant 1 .4 < z < 2.5 massive star forming galax- 
ies in the GOODS fields, described in a companion paper 
(Daddi et al. 2007, Paper I hereafter). As discussed in Paper I, 
we have built z ~ 2 galaxy samples free of detectable AGN 
activity. A sizable fraction shows distinct mid-IR excess with 
respect to that expected from star formation activity alone. 
This property of z ~ 2 galaxies implies that the extremely ef- 
ficient and penetrating 24 /im photometry from Spitzer is not 
always a reliable estimator of the ongoing star formation rates 
(SFRs) in galaxies, especially in the case of sources with very 
bright mid-IR luminosities (particularly those encountered at 
z ~ 2 in surveys shallower than GOODS). In this paper, we 
investigate the nature of these mid-IR excess sources in de- 
tail. Stacking in the X-ray bands with Chandra clearly un- 
veils heavily obscured AGNs, in many cases (or perhaps most 
cases) with Compton thick central BHs, which might be re- 
sponsible for the mid-IR excess. The 24 /im photometry is 
therefore still extremely valuable also in the case of bright 
24 /im galaxies at z = 2, as it allows us to investigate this 
previously unidentified population of distant AGNs, relevant 
for the history of accretion onto super-massive BHs and their 
feedback on galaxy formation and evolution at a crucial epoch 
of massive galaxy formation. 

The paper is organized as follows: in Section|2]we describe 
the mid-IR excess galaxy sample and summarize all relevant 
background work described in full detail in the companion 
paper (Paper I). Section [3] compares the mid-IR to UV es- 
timated SFRs (i.e., the mid-IR excess) to a number of other 
observed galaxy properties and shows that a warm dust com- 
ponent in addition to the colder star formation component 
is present in these galaxies at near-IR to mid-IR rest-frame 
wavelengths. The X-ray stacking analysis, unveiling the pres- 
ence of heavily obscured AGNs, are presented in Section |4j 
with the constraints on the absorbing gas column density sug- 
gesting that these sources are typically Compton thick. We 
discuss the implications of these results with regard to the ob- 
scured AGN fraction at high redshifts, the X-ray background, 
the coeval growth of BH and galaxies, and AGN duty cycle 
in Section |5l the possible relevance of these results for AGN 
feedback is discussed in Section [6] Future prospects for fur- 
ther developments in understanding this previously unknown 
population of z ~ 2 AGNs are given in Section [7] Summary 
and conclusions are provided in Section [8] We quote stellar 
masses and star formation rates for the case of a Salpeter IMF 
from 0. 1 and 100 M G , and we adopt a WMAP cosmology with 
fl A ,n M = 0.74,0.26, and h = H [km s" 1 Mpc^l/lOO = 0.73 
(Spergel et al. 2003). Unless explicitly stated otherwise, we 
quote magnitudes in the Vega scale. 

2. THE SAMPLE 

2.1. Galaxy selection and datasets 

The dataset and galaxy samples used in this paper are fully 
described in Paper I. For completeness we provide the rel- 
evant information here. We have selected sources in the K- 
band to completeness limits of K — 20.5 in the GOODS-North 
and K = 22 in the GOODS-South field, and then used the BzK 
color selection technique of Daddi et al. (2004b) to identify 
a total of ~ 1200 galaxies at 1.4 & z ~ 2.5. Accurate pho- 
tometric redshifts (a z ~ 0.25) as well as a large number of 
spectroscopic redshifts have been used. We exclude from the 
analysis the small fraction of contaminant BzK galaxies with 
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z < 1.2 or z > 3. Sources of spectroscopic redshifts are de- 
scribed in Paper I and include all available public datasets, 
plus currently unpublished redshifts from the ultradeep (15- 
30 hours per mask) spectroscopy of the GMASS survey (Kurk 
et al. 2007, in preparation). We have excluded AGN-hosting 
galaxies by eliminating from the sample (1) sources with a 
2-8 keV band detection in the 1 Ms and 2 Ms Chandra Deep 
Field X-ray catalogs of Alexander et al. (2003a) for GOODS- 
S and N respectively, and (2) sources with power-law SEDs 
over Spitzer IRAC and MIPS, lacking a distinct decrease in 
the flux density (/„) at wavelengths beyond 1 .6/im rest-frame, 
based on visual classification (eventually, this removes only 6 
galaxies in total due to power-law SEDs that were not also 
detected in the hard X-ray band). We also excluded 10% of 
the galaxies in the sample because of source blending due to 
the low resolution (~ 1 .6") at IRAC wavelengths, which gen- 
erally hampers measurement of accurate Spitzer MIPS 24 /im 
flux densities as well. 

The main aim of Paper I is to investigate the nature of star 
formation at high redshift and to identify the best mean to 
estimate the ongoing SFR activity. The work thus mostly fo- 
cuses on BzK star forming galaxies (as opposed to BzK pas- 
sive galaxies). All the multiwavelength star formation rate 
(SFR) indicators available in GOODS have been used. These 
include: (1) the deepest current Spitzer MIPS 24 /im imaging 
(Chary et al. 2007, in preparation), reaching about 12/iJy (3a) 
in both fields; (2) the deepest current Spitzer MIPS 70 /im 
imaging reaching about 1.8mJy in part of both GOODS fields 
(Frayer et al. 2006 and in preparation); (3) 160 hours of VLA 
1.4 GHz data in GOODS-N with 4.7/Jy of rms per beam 
(Morrison et al. 2007, in preparation); (4) 0.5-2 keV X-ray 
imaging from Chandra (used mainly in this paper, where 
we also discuss the AGN contamination that is needed for 
the interpretation of those data); (5) a deep 850 /im map of 
GOODS-N (Pope et al. 2005); and (6) UV photometry from 
GOODS ground-based and HST imaging (Giavalisco et al. 
2004). The way in which photometry at these various wave- 
lengths is interpreted to measure SFR in galaxies, and the im- 
plied known limitations, are described in Section 3 of Paper I. 
The main results of Paper I, relevant to this work, are pre- 
sented in Section 4 and 5 of that paper. In summary, we 
find that a fraction of the z = 2 galaxies (generally those with 
the brightest mid-IR luminosities) have significant excess flux 
density in the mid-IR at 24/Ltm, i.e. at rest-frame ~ 8/im, 
compared to the flux density expected on the basis of their 
ongoing star formation activity. A consistent mid-IR excess 
is also inferred from a comparison with radio, 70 /im, UV, 
submm (Paper I) (and soft X-ray, as discussed later in this pa- 
per) derived SFRs. We also find that the SFRs estimated from 
the dust reddening-corrected UV luminosities compares fairly 
well with respect to all other available indicators. 

The aim of this paper is therefore to carefully identify and 
study the population of sources with mid-IR excess. In or- 
der to work with meaningfully defined galaxy samples, being 
able to compute accurately SFRs from both mid-IR and from 
the UV, we limited in the following the analysis to 24 /im de- 
tected sources with well defined UV slopes (i.e. those with an 
error in (B — z) color less than 0.4 mag). This is required in or- 
der to obtain reasonable estimates of dust corrected UV SFR 
(SFRuy^on), with a formal error in the reddening correction 
lower than a factor of w 2. 

We further excluded galaxies that appear to be quiescent, or 
have only a low level of star-formation activity. These sources 
are generally very red in the UV because of relatively old stel- 



lar populations affecting the colors, and not simply because of 
dust reddening. For such sources the dust extinction correc- 
tion is expected to fail, leading to overestimated SFR\j\ con: , 
so that the comparison of mid-IR based SFRs to UV ones is 
not feasible. Based on the Paper I finding of a tight correlation 
between SFR and stellar mass at redshift z ~ 2 (as well known 
now to exist from z — to z — 1, Elbaz et al. 2007; Noeske et 
al. 2007), we excluded quiescent galaxies as those for which 
the specific star formation rate SSFR m j c i_iR < SSFR® njdIR /3, 

where SSFR = SFR^-ir/Moss and SSFR® ijd IR is the me- 
dian for 24 /im detected z ~ 2 galaxies in GOODS. This cri- 
terion also requires a 24 /im detection. 

2.2. Definition of mid-IR excess sources 

Following the results from Paper I, we define here mid-IR 
excess sources as sources satisfying 

Log(SFR mid _ IR+uv /SFR VVxon ) > 0.5 (1) 

Here Sf7? m!£ /-/fl+£/y is derived from the 24 /im flux density at 
the measured spectroscopic or photometric redshift, and us- 
ing the Chary & Elbaz (2001) templates, which are based on 
local mid- to far-IR correlations, with the addition of SFR 
estimated from the directly observed UV luminosity (non- 
reprocessed light from star formation). The term SFR\j\ con: 
is estimated from the UV after correcting for dust extinction 
based on the Calzetti et al. (2000) law (the average correction 
is a factor of ~ 40 for K < 20 sources, Daddi et al. 2005b, 
and decrease to a factor of a few to the faintest K = 22 lim- 
its). A galaxy can satisfy this mid-IR excess criterion either 
from overestimating the SFR m j c (_[R respect to SFRyy, or from 
underestimating the latter relative to the former. The first 
case corresponds to a genuine excess flux density in the mid- 
IR; the second one can be due to dust obscuration exceeding 
the Calzetti law prediction (optically-thick UV emission). In 
Paper I it is shown that the second case is relevant at most for 
a small fraction of the sample ( <, 15%), and therefore the vast 
majority (but not all) of mid-IR excess sources are overlumi- 
nous, for some reason, at mid-IR wavelengths. 

It is important to recall here that the luminosity dependent 
Chary & Elbaz (2001) templates that we adopt have been built 
on the basis of the local correlations in the mid-IR to far-IR 
for star forming galaxies. A mid-IR excess with respect to 
these SEDs is equivalent to a mid-IR excess with respect to 
the SED of standard star forming galaxies in the local uni- 
verse. In any case, it should be clarified that this population 
of mid-IR excess galaxies is well defined independently of the 
class of SEDs used to interpret the mid-IR luminosities. As 
shown extensively in Paper I, z = 2 galaxies at a given intrin- 
sic SFR or bolometric luminosity appear to display a wide 
range in their L(8/im). Therefore mid-IR excess galaxies can 
be also seen as galaxies with L(8/im) larger than the average 
< L(8/xm) > of galaxies of a given SFR (or of local star form- 
ing galaxies with the same SFR). Deviations from the normal 
L(8/im) versus SFR relation are also observed in the local uni- 
verse, and generally ascribed to the presence of an AGN (e.g., 
Genzel et al. 1998; Laurent et al. 2000; Dale et al. 2006; 
Armus et al. 2007). 

The threshold chosen in Eq. 1 (0.5 dex, about a factor 
of 3) is dictated by the evidence that, in the opposite di- 
rection to mid-IR excess galaxies, there are very few galax- 
ies with Sf7?uv,corr in excess of more than a factor of 3 of 
SFR m id-iR+uv (FigfU). This suggests that samples of gen- 
uinely discrepant galaxies with too high ratios are to be found 
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FIG. 1. — The ratio of mid-IR-based to UV-based star formation rate 
(SFR) is plotted as a function of the redshift. Here and in several of the 
figures in this paper, we plot only the GOODS-S portion of the galaxy sample, 
being deeper and extending to K < 22. Filled squares indicate spectroscopic 
redshifts, empty squares indicate photometric redshifts. The horizontal line 
defines the separation between mid-IR-excess and normal galaxies, as given 
by Eq. 1 . The rightmost panel shows the distribution of sources as a function 
of mid-IR excess ratio. The dotted line is the reflection of the bottom part 
of the histogram around a ratio of 1 . The error bar in the top-left part of the 
figure shows the typical error in the SFR ratio, inferred from the spread of the 
histogram around a ratio of 1 . 



above a similar threshold, assuming that intrinsic scatter in 
the ratio is to first order symmetrical. 1 We will later show that 
this threshold is well defined also with respect to the X-ray 
properties of the samples. 

Although the use of UV as the primary SFR indicator is less 
than optimal because of the possibility for substantial obscu- 
ration from thick lines of sight, it is the only one available for 
nearly all 24 /im sources at z ~ 2 in GOODS (the other SFR 
indicators being limited by sensitivity and applicable individ- 
ually to only a fraction of the galaxies, e.g., 20% for radio). 

The combination of the requirements discussed in the pre- 
vious section leaves us with a total sample of 430 (183) galax- 
ies in the GOODS-S (GOODS-N) field, or about 50% of 
the primary sample. Of these, 106 (58) satisfy Eq. 1 and 
are therefore classified as mid-IR excess galaxies — repre- 
senting 24.6% (31.7%) of the total sample. Likewise, robust 
spectroscopic redshifts are available for 108 (44) galaxies, or 
25. 1% (24.0%) of the total sample. Most of the mid-IR excess 
sources have very high 8 /mi rest-frame luminosities, often in 
excess of 10 Lq (FigO, and the vast majority of galaxies 
in our sample having L(8/im)> lO n L also shows a mid-IR 
excess. This is reminiscent of finding that most luminous IR 
galaxies in the local universe tend to be AGN dominated (e.g., 
Tranetal. 2001). 

3. MID-IR EXCESS AND OTHER GALAXY PROPERTIES 

1 It is worthwhile to emphasize that, despite the large dust reddening cor- 
rections to SFRuv.corr. most sources in Fig.[T]cluster around a ratio of 1, and 
the Gaussian component in the distribution has a dispersion of only ~ 0.2 dex 



With the aim of understanding the nature of the sources 
with mid-IR excess, and therefore the origin of this excess, 
we have searched for the existence of possible correlations of 
this excess with other galaxy properties. 

3.1. UV slopes, redshift, morphology and colors 

As shown in Fig. 13 of Paper I, the mid-IR excess does 
not strongly depend on the UV slope of the galaxies, and we 
therefore exclude that it is due to effects inherent in the cor- 
rection for reddening of the UV luminosity. A possible origin 
of the excess could instead be mid-IR SED evolution with 
redshift, within the 1 .2 < z < 3 range populated by our sam- 
ple. We have checked if mid-IR excess sources preferentially 
appear at some high or low redshift within the z ~ 2 sam- 
ple (Fig. [TJ. There is a tendency for z ~ 1.5 sources to have 
low mid-IR to UV ratios. This could be due to 9.7/im sil- 
icate absorption entering the 24 /im Spitzer MIPS bandpass, 
or to an overestimation of the reddening correction at these 
lowest redshifts. There is also a tendency for sources with the 
strongest excess to be located at z > 2. However, the median 
ratio of mid-IR to UV based SFRs is fairly stable as a func- 
tion of redshift, and no strong trend with redshift is detected. 
If we are seeing galaxies with anomalous mid-IR SEDs, com- 
pared to z ~ 0, they are rather homogeneously distributed in 
the probed redshift range. 

We have performed Sersic index fitting using Galfit (Peng 
et al. 2002) to the ACS F850W band images of our sample 
of z = 2 galaxies in GOODS-S. Results and full details will 
be presented in Ravindranath et al. (in preparation). Compar- 
ing the Sersic index and sizes distributions of mid-IR excess 
and normal galaxies, we do not find statistically significant 
differences for the full sample with K < 22 in GOODS-S. 

Instead, we find a strong correlation of the mid-IR excess 
with the (K — IRAC) colors. Fig. [3] shows that mid-IR ex- 
cess objects have redder (K — 5.8/im) colors. Galaxies with 
spectroscopic redshift follow the same trend of those with 
only a photometric redshift. Similar trends are observed for 
the other 3 IRAC channels (3.6 /im, 4.5 /im and 8.0 /im, not 
shown here). In principle, redder colors could suggest higher 
than average redshifts due to /^-correction effects, but this is 
ruled out by the weak trend with redshift that we observe, 
e.g., for the (K — 5.8/im) color. We conclude that the optical 
to near-IR rest-frame colors of mid-IR excess sources are in- 
trinsically redder than that of normal sources. These results 
imply that whatever the cause of 24 /im flux density excess 
in these galaxies, this excess is detectable also toward shorter 
wavelengths as seen in the IRAC bands. The distribution in 
the SFR ratio shown in Fig. [3] for K — 5.8/im< 1 colors sup- 
ports our choice of the factor of ~ 3 ratio in Eq. 1 for defining 
mid-IR excess sources. 

3.2. Median SEDs of galaxies with and without mid-IR 
excess 

The systematic difference observed for colors reflects dif- 
ferent typical SEDs between the UV/optical and the near- 
IR rest-frame. We have derived multiwavelength rest-frame 
SEDs of mid-IR excess and normal galaxies by computing 
the median magnitude for galaxies in each sample. Results 
are shown in Fig.g]for GOODS-N and GOODS-S. This visu- 
ally confirms that excess emission is detected at all rest-frame 
wavelengths greater than 1/im. The typical SFRu\, CO rr of the 
2 samples are very close for GOODS-N, as can be inferred 
from the similarity of the SEDs below 1/im [implying also 
similar median E(B — V)]. Being much deeper in /T-band, the 
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FIG. 2. — The mid-IR excess is observed to correlate with the mid-IR 
luminosity of galaxies (Paper I; the lack of sources in the upper-left part of 
the plot might be in part a selection effect). Symbols are like in Fig. [T] The 
majority of the sources with L(8/jm)> 10"L© would be classified as mid-IR 
excess sources based on the definition in Eq. 1 . Multiwavelength analysis in 
Paper I (and later in this paper) shows that the excess is present with respect 
to radio, UV, 70 (im, soft X-ray, and submm SFR indicators. 



GOODS-S mid-IR-normal sample has lower &F/?uv,corr than 
the mid-IR-excess sample by some 25-30%. Mid-IR excess 
sources tend to be found more often at brighter T^-mags, as 
one can argue from the figure; see also the next section. This 
explains the larger difference in the normalization of the SEDs 
above 1/im rest-frame. At the shortest wavelengths probed 
here, it is instead the normal galaxies that have brighter me- 
dian magnitudes. This is due to the fact that fainter galaxies 
are bluer (see e.g. also Daddi et al. 2004b). 

More importantly, we notice that the peak of the stellar 
dominated component of the SED (say, below 2-3/im rest- 
frame) of mid-IR excess objects is located at longer wave- 
lengths than in normal galaxies (e.g., IRAC 5.8/im band ver- 
sus 4.5/im band), where the peak is consistently located at 
~ 1.6/im as expected (e.g. Sawicki 2002). This difference 
persists when accounting for the slightly higher redshifts of 
mid-IR excess objects (z ~ 2.1 versus z ~ 1.9 on average). In 
order for this difference to vanish, one might assume that mid- 
IR excess objects are at even higher redshifts than inferred 
here. In turn, the mid-IR excess would vanish if we assume 
that these galaxies are actually systematically at lower red- 
shifts than what estimated here (which would increase more 
strongly SFR^-ir+uv than SFR\jy corr ). However there is 
no way to eliminate both the mid-IR rest-frame excess and 
near-IR rest-frame excess by simply moving the SEDs in red- 
shift. 2 This is an important piece of evidence demonstrating 

2 If these were otherwise normal galaxies (also in their mid-IR) with just 
photometric redshifts (and thus SFR mil j_i R+u v) largely overestimated, we 
should have found the peak of the stellar SEDs at lower rest-frame wave- 
lengths than in normal galaxies. If we had instead underestimated the red- 
shifts of the sources (but notice that objects with Zspec ar >d Zphot behave con- 
sistently in Fig. |5J, placing them at even higher redshifts would further in- 



FlG. 3. — The mid-IR to UV based SFR ratio is plotted as a function of 
the (AT - 5.8/jm) color, in the AB scale, for the GOODS-S field. Symbols are 
as in Fig. [T] The figure shows that mid-IR excess sources are systematically 
redder in K — 5.8/^m color than normal galaxies. 



that we are finding mid-IR excess galaxies not simply as the 
result of biases in the photometric redshifts. 

In principle, an SED peak longward of 1 .6/im could be due 
to extremely strong dust reddening (e.g., Chary et al. 2007). 
However, this is disfavored by the comparison of the UV 
slopes, and further by the agreement between the radio and 
UV SFRs which rules out substantial additional obscuration 
in the UV (beyond the obscuration that we estimate based on 
the UV slopes). We also tested that these results are robust if 
restricted to sources with spectroscopic redshifts only, and if 
excluding (in GOODS-N) sources that are radio detected and 
could be truly UV obscured cases. We conclude that from 
the SEDs there is compelling evidence that the mid-IR excess 
is intrinsic and extends to near-IR rest-frame wavelengths, in 
such a way that the SED shape is altered and its peak moved 
to wavelengths longer than 1.6/im (to about 1.8-1.9/im). The 
required excess component increases sharply with wavelength 
up to ~ 10/im rest-frame. 

Stacking of sources at 70 /im (after excluding the few indi- 
vidual detections) yields /t 3er detections for both mid-IR ex- 
cess and normal galaxies in both GOODS-N and GOODS-S 
fields (Fig. [4]). In both fields, mid-IR excess and normal galax- 
ies have 70 /mi stacked flux densities in agreement within the 
errors, although in both fields the mid-IR excess sources are 
marginally brighter (less than ler significance in both cases). 
This is consistent with the two classes having comparable un- 
derlying SFRs. This also implies that the excess emission 
is getting less important, in relative terms, when sampling 
~ 20/mi rest-frame wavelengths with the 70 /im imaging. The 
observed 70 /im to 24 /im ratio is 3.2 and 3 for mid-IR excess 
galaxies in GOODS-N and GOODS-S respectively, less than 
half the ratio expected for a local ULIRG if placed at z ~ 2.1. 
The ratios for the normal galaxies are in better agreement with 

crease the strength of the excess seen at 24 fim. 
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FIG. 4. — Median SEDs of mid-IR excess objects (black) compared to those of galaxies with more normal mid-IR properties (red). Left panel (a) is for the 
GOODS-N field, limited to K < 20.5; right panel (b) is for GOODS-S, to K < 22. The GOODS-S photometry includes BVIzJHK, the 4 Spitzer IRAC bands, 
and Spitzer MIPS 24 fim and 70 fim. In GOODS-N we have the same, plus U-band. 



those expected from the Chary & Elbaz (2001) models based 
on local correlations. This trend is similar to what found by 
Papovich et al. (2007). 

The mid-IR excess therefore manifests itself as a relatively 
warmer continuum dust emission, with respect to the star for- 
mation components, whose contribution peaks in the mid-IR 
but it is still detectable in the near-IR rest-frame. This kind 
of warm emission is unlikely to be star formation dominated, 
as the reprocessed emission from star formation declines very 
steeply below 5fjm (Genzel et al 1998; Laurent et al 2000), 
and is a clear signature of the presence of an AGN. These 
galaxies could be similar to those selected as putative AGN 
hosts on the basis of power-law continua over the same wave- 
length range (e.g., Alonso Herrero et al. 2006; Donley et 
al. 2007), but with a smaller relative contribution of the AGN 
at least in the IRAC bands, at near-IR rest-frame wavelengths. 
By comparing the excess at near-IR and mid-IR rest-frame 
wavelengths, we find that the additional contribution has a 
shape comprised between that of the nucleus of NGC 1068 
and of Mrk 231 (at a fixed 24 /zm excess, NGC 1068 under- 
predicts and Mrk 23 1 overpredicts the median excess over the 
IRAC bands). 

This contamination observed in the IRAC bands will have 
consequences for SED fitting with stellar population models 
to derive stellar population parameters and stellar masses. We 
explore this issue in more detail in Maraston et al. (in prepa- 
ration). This contamination can also affect the estimation 
of photometric redshifts of z ~ 2 galaxies when using IRAC 
bands, as this estimate relies also on the location of the 1. 6 fim 
bump of the stellar SEDs. Comparing to our set of spectro- 
scopic redshifts, we can in fact detect a small bias, with the 
photometric redshifts of mid-IR excess galaxies being slightly 
overestimated. Compared to normal galaxies, the median dif- 
ference between photometric and spectroscopic redshifts is 
larger by about 0.04 (0.09) in GOODS-S (GOODS-N), just 
a few percent of ( 1 + z) ■ 



It is also important to recall here that the warm continuum 
contamination of the SEDs disappear, or at least is entirely 
negligible, over the UV rest frame. This is based on the me- 
dian SEDs shown, and also on the consistency of UV based 
SFR estimates and submm, radio, and 70 /im, for both mid- 
IR excess and normal galaxies, as extensively discussed in 
Paper I. In addition, no AGN emission lines are generally de- 
tectable in the UV from the mid-IR excess objects, even in the 
ultradeep 15-30 hour integration GMASS spectroscopy. This 
is also consistent with the result that the rest-frame UV mor- 
phological properties of mid-IR excess and normal galaxies 
are similar. 

3.3. Mid-IR excess and stellar mass 

We have checked if and how the incidence of mid-IR excess 
galaxies varies with the inferred stellar mass (Fig. [5j. Stel- 
lar masses for our sample are estimated from the empirically 
calibrated relations presented in Daddi et al. (2004b), as dis- 
cussed in Paper I, based on the observed photometry in the B-, 
z- and K- bands. We find that mid-IR excess sources tend to 
preferentially reside in galaxies with larger than average stel- 
lar masses. The fraction of mid-IR excess sources increases 
steeply with the mass of the galaxies. 

We have investigated if the relation shown in Fig.|5]could be 
the result of a selection effect. First, we have checked that the 
trend formally remains if we add to the sample the galaxies 
with uncertain measurements of the UV slope (and therefore 
of the UV SFR and mid-IR excess). Also, the sample plotted 
in Fig. [5] is not mass complete to the lowest masses it for- 
mally reaches (~ 0.5 x 10 Mr.,). Given our K-band selection 
to K < 22 (for the GOODS-South sample shown in the fig- 
ure on which this analysis is based), we expect to be largely 
complete only above ~ (2-3) x 1O 1O M . We would miss star 
forming galaxies below this stellar mass threshold if they have 
very red colors. However, we have verified that we still find 
statistical evidence for an increasing fraction of mid-IR excess 
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FIG. 5. — Left panel (a): the ratio of mid-IR versus UV inferred SFR is plotted against stellar mass. Colors and symbols are as in Fig. fj] Right panel (b): the 
fraction of mid-IR excess objects inside star forming galaxies as a function of stellar mass. Bins are spaced by 0.4 dex and contain sources within ±0.2 dex in 
stellar mass (i.e., the 5 bins plotted are independent). The line is a weighted linear lit to the data. Error bars on the fractions are Poisson, while error bars in the 
x-axis reflect the width of the bin. 



sources with stellar mass if we limit the sample to the range 
in stellar masses that is expected to be highly complete. Sim- 
ilarly, if we limit our analysis to blue galaxies only, for which 
we are mass complete over nearly the whole mass range, we 
still recover the trend shown in Fig. [5] 

The warm continuum contribution detected over the IRAC 
bands is unlikely to bias the results on the dependence of the 
fraction of mid-IR excess galaxies with stellar mass. Already 
over the IRAC bands the typical excess contribution is less 
than 0.5 mags 3 and we recall that the stellar mass estimates 
are not derived using the IRAC bands, but only up to TT-band. 
The presence of a few tenths of mags of AGN contribution in 
the observed Zf-band would not alter substantially the trend. 
Therefore we conclude that this result is solid and not simply 
a selection bias. 

In Paper I, we describe a strong correlation between stellar 
mass and SFR\]\ colT . Because of this correlation, the fraction 
of mid-IR excess objects also increases with SFR\jy COTr . This 
effect is less pronounced and less evident than the one found 
for the stellar mass. 

4. X-RAY PROPERTIES 

The overall properties of mid-IR excess sources, as ex- 
plored in the previous sections, point to the presence of a 
warm mid-IR continuum component, which could be due to 
AGN activity. Since we have removed from the sample all 
X-ray identified AGNs in the 2-8 keV band, any powerful 
(Lx > 10 43 erg s _1 ) AGN present in the mid-IR excess ob- 
jects would have to be heavily obscured in order to escape 
being individually detected in the X-rays (note that the ob- 
served 2-8 keV band corresponds to the rest-frame 6-24 keV 

3 This is true for both GOODS fields; we recall that the larger difference in 
the SED normalizations over the IRAC bands that is observed in GOODS-S 
is mainly due to the different median stellar mass between mid-IR excess and 
normal galaxies (Fig.[5J. 



band). We can search for the signature of heavily obscured 
AGNs below the X-ray detection limit using X-ray stacking 
analyses (e.g., Brandt et al. 2001). Here we use the code 
adopted by Worsley et al. (2005), which stacks sources and 
calculates the significance of the stacked result using 10,000 
Monte Carlo trials. In all cases, we limit the X-ray stacking to 
sources within 5.5 arcmin of the Chandra aim point, to maxi- 
mize sensitivity. 

4. 1 . X-ray stacking results 

We stacked 59 mid-IR excess galaxies in GOODS-S and 
found a 5a detection in the 2-8 keV band and an 8er detec- 
tion in the 0.5-2 keV band (Fig. |6|. The count ratio (H/S) 
is 0.90 ±0.13, corresponding to a nominal power law pho- 
ton index T = O.8+Q3. Using T = 0.8, we infer an average 

2-8 keV flux of 5.6 x 10~ 17 erg s _1 cm~ 2 and an average 
0.5-2 keV flux of 1.1 x 10~ 17 erg s" 1 cm" 2 . The flat X-ray 
spectral slope strongly argues for the presence of heavily ob- 
scured AGN (e.g., see Fig. 2 of Alexander et al. 2005b). The 
175 normal GOODS-S sources (i.e., those without mid-IR ex- 
cess) give a stacked 3er signal in the 2-8 keV band and a 12<r 
detection in the 0.5-2 keV band (Fig. [6]). The band ratio is 
0.36, corresponding to a spectral index r ~ 1 .7, which is fully 
consistent with star formation dominating the X-ray emission 
(e.g., Persic & Raphaeli 2002). Using this value of T we de- 
rive a 2-8 keV flux of 1.6 x 10~ 17 erg s cm~ 2 and a 0.5- 
2 keV flux of 1.0 x 10~ 17 erg s cm~ 2 . Stacking analysis 
of the smaller GOODS-N sample, that is limited to K < 20.5, 
gives a slightly less secure detection and lower S/N ratios, al- 
though still with hard X-ray fluxes twice as large for mid-IR 
excess versus normal galaxies. Similar properties are seen in 
the GOODS-S sample when limited to K < 20.5. In the fol- 
lowing we concentrate on the GOODS-S sample to K < 22, 
because it has the highest S/N X-ray detections. The results 
of the X-ray stacking are listed in Table 1 . 
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TABLE 1 . X-RAY STACKING ANALYSIS OF GOODS z = 2 GALAXIES 



SAMPLE 


N 


Exp 


Counts 


Counts 


Counts 


S/N 


S/N 


S/N 


Flux 


Flux 


Flux 






SB 


FB 


SB 


HB 


FB 


SB 


HB 


FB 


SB 


HB 






Ms 


10~ 6 s _I 


10~ 6 s _I 


10 6 s 1 








10 17 cgs 


t n 17 

10 17 cgs 


10" 17 cgs 


NORMAL: 
























K < 20.5 N 


69 


121.2 


2.05 


1.42 


0.65 


6.6 


8.4 


2.4 


1.78 


0.74 


1.32 


K < 20.5 S 


58 


49.3 


3.33 


2.35 


1.01 


6.8 


8.4 


2.5 


2.94 


1.25 


2.24 


K < 22.0 S 


175 


146.7 


2.52 


1.86 


0.68 


9.0 


11.7 


2.9 


2.23 


0.99 


1.50 


EXCESS: 
























K < 20.5 N 


27 


47.2 


3.80 


2.57 


1.26 


7.3 


9.0 


2.8 


3.29 


1.35 


2.55 


K< 20.5 S 


27 


22.8 


4.09 


2.88 


1.24 


6.0 


7.5 


2.2 


3.61 


1.53 


2.76 


K < 22.0 S 


59 


49.8 


4.01 


2.13 


1.92 


8.6 


8.0 


5.0 


6.84 


1.07 


5.62* 



NOTE. — N is the number of sources that were stacked. FB is full band (0.5-8 keV), SB is soft band (0.5-2 keV), and HB is hard band (2-8 keV). 10000 Monte-Carlo trials are 
performed to estimate the background and S/N for each sample. Effective exposure time is measured in the soft band. Fluxes are calculated using the X-ray spectral slope of T = 0.8 
for the K < 22 GOODS-S mid-IR excess galaxies and assuming T — 1.9 for the samples that are not detected in the hard band. 
* using the model shown in Fig.[9]we derive a ~ 1.7 X higher flux of 9.55 x 10~ 17 erg s -1 . 





FIG. 6. — The soft and hard X-ray Chandra stacked images of normal 
and mid-IR excess z = 2 objects in GOODS-South. Left panels show the 
soft (0.5-2 keV) bands, while the right panels are for hard (2-8 keV) bands. 
The top 2 images are for normal galaxies, while the bottom 2 are mid-IR 
excess objects. It is evident that similar soft X-ray fluxes are detected in the 
2 samples, but much stronger hard emission is detected for the mid-IR excess 
z ~ 2 galaxies. Images have been smoothed by a Gaussian with the size of 
the PSF. The circles (4" diameter) show the expected location of the signal. 



The observed X-ray fluxes can be directly used to esti- 
mate luminosities, without accounting for obscuration. For 
the mid-IR normal galaxies in the GOODS-S sample (< z >= 
1.9), we infer isotropic luminosities of 2.7 x 10 erg s _1 in 
the 2-8 keV rest-frame band (mapping very closely to the 0.5- 
2 keV observed band) and of 4.2 x 10 41 erg s in the 5.8- 
23.2 keV band rest-frame, assuming a T = 1 .9 spectral slope. 
For mid-IR excess galaxies in GOODS-S (< z >— 2.1) the 
fluxes correspond to X-ray luminosities of 3.8 x 10 41 erg s _1 
in the 2-8 keV rest-frame band and 2 x 10 42 erg s _1 in 
the 6.2-24.8 keV band rest-frame, for the spectral slope of 
r = 0.8 derived from the count ratios. 

We have explored how these X-ray stacking results 
depend on the threshold of 0.5 dex in the ratio of 
SFR(mid - IR + UV) to SFR(\JV, corr) that we have adopted 
in Eq. 1 for defining mid-IR excess sources. We have varied 
the threshold from 0.2 dex to 1.0 dex, and performed X-ray 
stacking separately for sources above and below the threshold. 
Fig.|7]shows that the general X-ray properties of the "normal" 



populations are quite insensitive to the threshold at which the 
excess sources are identified, given that the latter ones are in 
all cases a minority. For the mid-IR excess sources, the X-ray 
stacking confirms that a threshold 2-5 in the SFR ratio effec- 
tively singles out the population of hard X-ray emitters (with 
similar count rates in both hard and soft Chandra bands). The 
X-ray count rates, and therefore the average luminosities of 
the mid-IR excess sources, increase with the threshold, while 
their space density obviously decreases, but the accuracy in 
the measurements declines due to lower-number statistics. 

A possible concern with stacking analyses is whether a 
small subset of bright sources dominate the stacked results. 
We verified that this is not the case by performing stacking in 
galaxy subsamples. The results shown in Fig.|7]argue that our 
X-ray stacking conclusions are robust to the exact definition 
of the sample. 

Although the X-ray band ratio indicates the presence of ob- 
scured AGNs, it does not fully characterize the properties of 
the sources. We can gain additional insight into the X-ray 
spectral properties of mid-IR excess sources by stacking the 
data inside narrower energy bands. Stacking the GOODS-S 
data using the sub bands 0.5-2 keV, 1-2 keV, 2-4 keV, 4- 
6 keV, 4-8 keV, and 6-8 keV, one obtains the results shown 
in Figs. [8] and [9] and Table 2. Significant ( ^ 3ct) detections 
are found in all of the sub bands, except for the 6-8 keV band 
where the detection is only at the 1.8ct level (still noteworthy 
given the low sensitivity of Chandra at these high energies). 
In summary, the mid-IR excess sources exhibit a hard X-ray 
spectrum which becomes harder with increasing energy, an 
unequivocal signature of highly obscured AGN activity. 

4.2. X-ray inferred SFRs 

The soft X-ray (0.5-2 keV band) fluxes of the samples 
(mapping quite closely to the rest-frame 2-8 keV luminosities 
at z = 2) can be used to estimate SFRs, since for star form- 
ing galaxies this is proportional to the population of massive 
X-ray binaries (Ranalli et al 2003; Persic et al 2004; Horn- 
schemeier et al 2005). The rest-frame 2-8 keV luminosity of 
mid-IR excess sources is ~ 30% higher than that of mid-IR 
normal sources, and would thus suggest ~ 30% higher SFRs 
for the formers. This is in excellent agreement with the es- 
timate obtained from the UV, where we find the same offset 
between the average SFRyy corr of the two samples and in the 
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TABLE 2. STACKED X-RAY SPECTRA 



1 10 
Threshold in SFR(mid-IR + UV)/SFR(UV, corr) 

FIG. 7. — (Top panel): X-ray counts from stacking of mid- 
IR excess versus mid-IR normal sources for a variety of thresholds 
separating these two classes: from 0.2 to 1.0 dex in the ratio of 
SFi?(mid-IR + UV)/SFfi(UV,corr). Black lines are for mid-IR excess 
sources (with ratio above the threshold), while red lines are for mid-IR nor- 
mal sources (with ratio below the threshold). Solid lines show hard (2-8 keV 
band) counts, while dotted lines show soft (0.5-2 keV band) counts, from 
stacking Chandra data in the lMs GOODS-S field. (Bottom panel): the sky 
density of the two populations, as a function of the defining threshold. We 
are not applying here completeness corrections to the sample, see Sect. 5.1. 



same direction. Interestingly, the proportionality that we de- 
rive is 



SFR m>con [M Q yr- 1 ] = 2.25 x 1(T 4U L 2 _ 8 keV [erg s~ l ] (2) 

which is very close to the relation given in Ranalli et al. (2003) 
(the Persic et al. 2004 relation only accounts for the contri- 
bution of high mass X-ray binaries). The agreement between 
the UV and X-ray estimated SFRs, for both samples of nor- 
mal and mid-IR excess galaxies in GOODS-S, suggests that 
no substantial AGN contribution is detected in the observed 
0.5-2 keV band of these mid-IR excess objects. Similarly, we 
find that the soft X-ray fluxes in the mid-IR excess samples in 
both fields, when limited to K < 20.5, are fully accounted for 
by the SFR inferred from UV. 

4.3. Evidence for Compton thick AGNs 

Estimating the obscuring column density for the AGNs in- 
side mid-IR excess sources is crucial for deriving their in- 
trinsic properties. Their stacked Chandra 0.5-2 to 2-8 keV 
band ratio implies an absorbing gas column density TVh ~ 
10 23 3 cm~ 2 , for an underlying standard AGN power-law con- 
tinuum with r = 1.9. Therefore, these AGNs are clearly ob- 
scured, but this estimate of Nn should be regarded as a strict 
lower limit because star-formation is contributing to the emis- 
sion, especially in the soft band. The 0.5-2 keV emission of 
the mid-IR excess galaxies can be entirely ascribed to their 
ongoing star formation, as discussed in the previous section. 
If we consider that the AGN component in the 0.5-2 keV band 
is consistent with zero given the errors, we formally obtain a 



BAND 


Counts 


S/N 


Flux 




10~ 6 s" 1 




10" 17 cgs 


0.5-lkeV 


0.57 


2.9 


0.34 


l-2keV 


1.67 


8.5 


0.79 


2-4keV 


0.89 


3.8 


1.42 


4-8keV 


1.16 


3.4 


5.25 


6-8keV 


0.47 


1.8 


4.94 



NOTE. — Results of the X-ray spectral stacking for mid-IR excess objects in 
GOODS-S to K < 22. Fluxes calculated using the overall X-ray spectral slope of 
r — 0.8, estimated from the X-ray band ratio. S/N of the data calculated using the 
background estimated from 10,000 MC trials. 



2(7 lower limit N\\ > 10 239 cm~ 2 for an intrinsic AGN power 
law continuum with T = 1 .9. 

We further used the result of spectral stacking shown in Ta- 
ble 2 and Fig. [9] to improve our constraints on 7%. The right 
panel of Fig. [9] shows the X-ray spectrum at 2-20 keV rest- 
frame energies of the mid-IR excess sources, after subtraction 
of a star-formation component (dashed line model shown in 
left panel of Fig. |9]). We have fitted this AGN component 
with the Gilli, Comastri & Hasinger (2007) models for ob- 
scured AGNs, where 7% is increased in steps of 1 dex up 
to Na = 10 25 5 cm~ 2 (this approach is admittedly crude, but 
there is little point for a more refined approach in this work). 
At the highest levels of obscuration, beyond the Compton 
thick regime of Nn = 10 24 cm~ 2 , reflection components start 
to be important in the AGN emission in the rest-frame ~2- 
8 keV (observed Chandra 0.5-2 keV for z ~ 2), as well as the 
emission from the Fe Ka line at 6.4 keV, and both are prop- 
erly included in the models of Gilli et al. (2007). Having con- 
volved the templates with the observed redshift distribution 
of our sample, the best fit is obtained for 7% = 10 2 



,24.5 cm -2 



cm 

,2 



and 7V H = 10 25 5 cm 



with both cases 7V H = 10 23 5 

producing a worse fit with A\ z > 5.8. According to Avni 
et al. (1976), this implies that they can be rejected at the 
> 98.5% confidence level. For the best fitting case of TYh = 
10 24 5 cm" 2 , the observed 2-8 keV band would be obscured 
by a factor of ~ 4, with respect to the intrinsic, unobscured 
emission. 

Therefore, the X-ray spectrum of mid-IR excess sources is 
fully consistent with the presence of two major components: 
a star-formation component (modeled here with a power-law 
emission of T = 1 .9), plus a highly obscured AGN component 
with (average) Na 



10 24 5 cm 2 , i.e., mildly Compton thick. 



4.4. Column density and mid-IR luminosities 

The column density of obscuring gas can be independently 
estimated by exploiting the correlations between mid-IR and 
X-ray luminosities for AGNs, in the assumption that the mid- 
IR excess is entirely due to the AGN. The mid-IR excess 
sources to K = 22 in GOODS-S have average (median) 24 /im 
flux densities of 145 (131) /Jy, while normal sources have 
average (median) 24 /im flux densities of 41 (27) /iJy, consis- 
tent with the ongoing SFR on the basis of local correlations. 
We have computed for all galaxies the 24 /im flux densities 
expected for their Sf7?(UV,corr) on the basis of Chary & El- 
baz (2001) models, and compared to the observed 24 /im flux 
densities. We find also in this way that for mid-IR excess 
galaxies the average (median) excess flux density is 100/iJy 
(80/iJy). The ratio of 24 /im excess flux density to the to- 
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FIG. 8. — Chandra sub-band stacking detections of mid-IR excess z ' 
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FIG. 9. — Left Panel (a): X-ray measurements of flux per unit energy (F(E)) for mid-IR excess sources (black) and normal sources (red). It is evident 
that mid-IR excess sources are much harder, inconsistent with pure star-formation as indicated by the photon index T = 1.9 power-law (Starburst). Sub-band 
stacking of mid-IR excess sources (black points with solid error bars) reveals that their F(E) continue to rise at the highest energies. These sources are consistent 
with a combination of star-formation plus a heavily obscured AGN component rising toward the highest energies. Right Panel (b): the result of subtracting the 
star-formation rate component from the sub-band Chandra stacking of the mid-IR excess galaxies, normalized to the soft band point (allowing for a plausible 
range of uncertainty to this quantity does not substantially alter the results). The overplotted line is an AGN model with N H = 10 24 5 cm~ 2 (Gilli et al. 2007) 
convolved with the redshift distribution of the sample. The bump in the models near 2 keV (in the observed frame) is due to the strong Fe emission line expected 
to be prominent in obscured AGNs. 



tal 24 /im flux density is found to be 69% (68%) on aver- 
age (median). Therefore, we infer that the typical excess 
flux density is w 80-100 /dy and that the AGN is contribut- 
ing typically some 2/3 of the observed 24 /im flux density in 
these mid-IR excess galaxies. These figures imply a ratio of 
24 /im flux density to hard X-ray flux of ~ 2 x 10 15 in the 
[mJy]/[erg s cm~ 2 ] units used by Rigby et al. (2004). Ac- 
cording to Rigby et al. (2004), based on the local AGN tem- 
plates from Silva, Maiolino & Granato (2004) with a range of 
obscuration, such large ratios at z = 2 would imply column 
densities A/h >> 10 cm~ 2 , or substantially Compton thick 
sources. Using the observed correlations between the AGN 
mid-IR light and unobscured hard X-ray light from Lutz et 
al. (2004), which appear to hold independently of obscura- 
tion at least to A/h ~ 10 25 cm~ 2 , we estimate that unobscured 
AGNs would contribute only « 3-10/iJy at 24 /im for z = 2 
given their 2-8 keV X-ray fluxes. In order to match the Lutz 
et al. (2004) correlation, the 2-8 keV observed flux requires 
a correction factor of ks 10-15, i.e., higher than the factor 
~ 4 derived in the previous section, and potentially implying 



Compton-thick AGN with A^h > 10 25 cirr 2 . 

It is instructive also to compare to the properties of 
NGC 1068, an extremely obscured local starburst harboring 
a heavily Compton thick AGN with Nh > 10 25 cirT 2 . If 
placed at z — 2, the NGC 1068 nucleus would have a hard 
band X-ray flux of 7 x 10~ 18 erg s _1 cm~ 2 and a 24 /im flux 
density of 25/iJy. This would give a mid-IR to X-ray ra- 
tio of ~ 3.5 x 10 15 ([mJy]/[erg s _1 cm -2 ]) fully comparable 
(slightly larger) than that of the mid-IR excess z — 2 galaxies. 
Therefore, these z = 2 sources are consistent with being the 
high-redshift analogs of NGC 1068, with overall luminosity 
scaled up by a factor ~ 5. The AGN emission of NGC 1068 
is thought to be obscured by over a factor of 10, even in the 
6-24 keV band (rest frame). 

We conclude from the mid-IR to X-ray comparisons that 
the majority of AGNs shrouded inside massive mid-IR ex- 
cess z = 2 galaxies are likely Compton thick, with average 
A/h > 10 24 cm~ 2 , perhaps even in excess of the estimate of 
A/h ~ 10 24 5 cm~ 2 derived solely from the X-ray stacking, and 
possibly in some cases they have extreme obscuration with 
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N H > 10 25 cm" 2 . 

4.5. Intrinsic AGN luminosities 

A determination of the intrinsic X-ray luminosity is crit- 
ical to assess the general implications of our findings, but 
unfortunately it is made difficult because the obscuration ap- 
pears to be very high. For the best fitting Gilli et al. (2007) 
mildly Compton thick model having A^h = 10 245 cm~ 2 , the 
observed 2-8 keV band is obscured by a factor of ~ 4. Con- 
sistently using this model to derive 2-8 keV flux from the 
observed counts yields a flux larger by a factor of 1.7 than 
that reported in Table 1, computed for T = 0.8. Correspond- 
ingly, we estimate that the average, isotropic, unobscured 
2-8 keV rest-frame luminosity of mid-IR excess galaxies is 
^2-8keV ~ 10 43 erg s -1 . To give an idea of the uncertainty in 
the correction for obscuration, we find that for the Gilli et al. 
(2007) templates the luminosities would decrease by a factor 
of 2.9 for Na = 10 23 5 cm~ 2 and would increase by a factor of 
3.5 for A^h = 10 25 5 cm~ 2 (although we emphasize that the ab- 
sorption correction for the case of A^h = 10 25 5 cm~ 2 is highly 
dependent on the assumed reflection efficiency). 

As a comparison, we can use the Lutz et al. (2004) corre- 
lation between unobscured X-ray luminosities and 6/im lu- 
minosities for AGNs, which is found to be independent of 
obscuration and has a scatter for individual sources of about 
a factor of 3. The relationship between the X-ray and mid- 
IR emission is strongly dependent on the dust-covering fac- 
tor around the central source and hence a source with a large 
dust covering factor will be more mid-IR luminous than one 
with a smaller dust covering factor (see Alexander et al. 
2005b). The average 24 /im excess flux density attributed to 
the AGN is converted to an average 8/im luminosity of order 
4 x 10 10 Z<g for the ensemble of mid-IR excess objects. We 
estimate ~ 70% this value for the typical 6/im luminosity of 
our sources (based on typical SEDs, see also Fig. 0). From 
Fig. 7 in Lutz et al. (2004) we would then estimate L2-8keV ~ 
4 x 10 43 erg s _1 , a factor of 4 larger than that obtained di- 
rectly from X-ray. In summary, plausible values for the typ- 
ical rest-frame 2-8 keV unabsorbed luminosities of mid-IR 
excess galaxies appear to be within (1-4) x 10 43 erg s _1 . 

This result can be cross-checked by looking at the X-ray 
luminosities of individually hard X-ray detected AGNs in our 
GOODS samples of BzK selected sources with 1 .4 < z < 2.5. 
If we are finding the heavily obscured, Compton thick, coun- 
terparts of the already known population of z ~ 2 relatively 
unobscured sources, we might expect that the Compton thick 
sources should have typical luminosities similar to the known 
AGNs, based on AGN unification models (Antonucci 1993; 
Urry & Padovani 1995). In fact, we find that the average 2- 
8 keV rest-frame luminosity of individually hard X-ray de- 
tected AGNs in our /f-selected sample is about Z,2-8keV ~ 
3 x 10 43 erg s _1 , neglecting obscuration corrections that are 
likely negligible in most of these cases. The individually X- 
ray detected z — 2 AGNs have a sky (and thus space) density 
that is a factor of ~ 3 smaller than that of our mid-IR excess 
galaxies in GOODS-S. Given the estimated X-ray luminosi- 
ties, the mid-IR excess galaxies are also consistent with the 
correlation between optical/near-IR and 2-10 keV rest-frame 
luminosities found by Brusa et al. (2005) for distant AGNs. 

One might also hope to constrain the bolometric luminosi- 
ties of the AGNs in our sample, although uncertainties in the 
AGN bolometric correction factors are currently large (see, 
e.g., discussions in Pozzi et al. 2007 and Hopkins et al. 2007). 



Assuming an AGN bolometric correction derived from the 
Elvis et al. (1994) SED, from the 2-8 keV X-ray luminosity 
we estimate Z^ol ~ ( 3 " 12 ) x l ° U er § s ~'' or about (°- 8 " 
3) x 1O U L0. The large range reflects only the uncertainty in 
the derived 2-8 keV rest-frame luminosity discussed above, 
as we are neglecting the unknown uncertainty in the bolomet- 
ric correction for these galaxies. If instead we start from the 
mid-IR excess, with a typical luminosity of 10 Lq, the Elvis 
et al. (1994) SED would predict a bolometric correction of 
a factor of 10, or ~ 1O 12 L . Given that our galaxies 
are heavily obscured in X-rays to a similar degree, it may be 
more appropriate to use the mid-IR bolometric correction for 
the NGC 1068 nucleus, which is a factor of 3. This would 
lead us to estimate Z^ol ~ 3 x lO 11 ^©- So, we are left with 
a plausible range spanning from 0.8 to 3 x 10 Z,©, although 
a value as high as 1O 12 Z,0 cannot be completely excluded at 
this stage. Admittedly, a very large uncertainty affects the es- 
timated AGN bolometric luminosity. 

From the ongoing star formation activity we estimate 
L S g OL ~5x 10 Lq for the typical mid-IR excess object, com- 
parable to the bolometric luminosity of the AGN. 

5. IMPLICATIONS 

5.1. The sky and space density 

The main result of our study is the uncovering of a large 
population of Compton thick AGNs at z ~ 2. We characterize 
here the abundance of this population. The sky density of the 
mid-IR excess AGNs in the GOODS-S field is ~ 0.7 arcmin -2 
to K < 22 (59 sources in a 5.5' radius region), after correcting 
for the 10% incompleteness in our sample due to the require- 
ment of non blended IRAC photometry (see Paper I). This 
estimate could increase by a factor of ~2, accounting for the 
sources that we excluded from the analysis not being able to 
properly measure for those either SFR(mid — IR + UV) or 
SFR(W,corr) in Eq. 1. (which requires a well measured 
UV slope and 24 /im detection, mainly, see Sect. 2). As the 
excluded sources lacking 24 /im detections are less likely to 
be mid-IR excess galaxies, we conservatively correct the sky 
densities only by a factor 1.3, to account for the overall frac- 
tion of red star forming galaxies that were excluded from the 
analysis for not having a well measured UV slope. This im- 
plies a sky density of 0.9 arcmin~ 2 , or some 3200 deg- 2 . Us- 
ing the volume within 1.4 < z < 2.5 (see Fig. 2 in Paper I) 
we infer a space density of 2.6 x 10 -4 Mpc~ 3 for the adopted 
ACDM cosmology. 

We are dealing with detections of average properties only 
(X-ray stacking), and there is likely a spread in luminosities 
and A^h among the sample. In principle only a fraction of 
the stacked mid-IR excess sources might be contributing the 
detected hard X-ray flux, e.g. in the case of a large con- 
tamination of the mid-IR excess galaxy sample from galax- 
ies with peculiar mid-IR SEDs that are otherwise powered by 
star formation. In that case, the space density of this obscured 
AGN population would be lower by the same fraction, and its 
average X-ray flux and luminosity would be proportionally 
higher. The contributing fraction cannot be too small, i.e. at 
least « 0.5, otherwise we would have detected directly most 
of these sources in our Chandra data. Also, the similarity 
in inferred X-ray luminosities of these z = 2 Compton thick 
AGNs to individually detected z — 2 sources (Sect. 4.5) argues 
against a lower fraction since this would boost their intrinsic 
luminosities. With this caveat, we have included all mid-IR 
sources in the computation of the space density of Comp- 



12 



E. Daddi et al. 




Redshift 

FIG. 10. — The space density of Compton thick AGNs. The blue circle is 
our z = 2.\ estimate, where we allow for a factor of 2 uncertainty. The red 
points are taken from Tozzi et al. (2006), based on direct Chandra detections 
in GOODS-S. The green points show the density that we crudely estimate 
for the survey of Polletta et al. (2006) and Martinez- Sansigre et al. (2006), 
accounting in the latter case for completeness correction due to their radio 
preselection. The lines show the predictions of the background synthesis 
model of Gilli et al (2007), as a function of the limiting X-ray luminosity. 



ton thick AGN candidates assuming that the mid-IR excess 
is due to the AGN as discussed above. Relaxing this hypothe- 
sis would not substantially affect any of the implications that 
we derive from this work. 

5.2. The predominance of Compton thick AGNs at high 
redshift 

It is interesting to compare the properties of this previously 
undetected population of obscured, Compton thick, AGNs 
to AGN populations already known. Over the same red- 
shift range 1.4 ,$ z 2.5, the sky density of mid-IR excess, 
Compton thick AGNs is a factor of ~ 2 higher than that of 
X-ray sources individually detected in X-rays in GOODS-S 
(Zheng et al 2004), a factor « 10 higher than X-ray detected 
submillimeter-emitting galaxies (Alexander et al. 2003b), and 
« 20 times higher than X-ray undetected power-law AGN 
candidates (Donley et al. 2007). 

We can assess the obscured-to-unobscured AGN ratio by 
comparing to the space density of unobscured AGNs se- 
lected in the soft band by Hasinger et al. (2005). For Lx > 
10 42 erg s -1 , they find a space density of ~ (3—4) x 10~ 5 at 
Z = 2, a factor of 5 smaller than what we find for obscured 
mid-IR excess AGNs at the same redshift. This implies that 
the ratio of obscured to unobscured AGNs is larger than at 
least 5, and is probably much larger yet, consistently with the 
models expectations (e.g., Gilli et al 2007). Comparing to the 
luminosity function for hard-X-ray selected AGNs at z = 2 
(Ueda et al. 2003; La Franca et al. 2005), the Compton thick 
sources that we have found have a space density that is larger 
by a factor of 2-3 than that of the sources already known and 
overall with similar luminosities. The difference would de- 



crease if it is only a fraction of mid-IR excess sources that are 
hosting an X-ray luminous AGN. A more detailed comparison 
is hampered by the uncertainties inherent the detailed prop- 
erties of our samples, as discussed in the previous sections. 
Overall, these results suggest that a dominant fraction of X- 
ray luminous AGNs are heavily obscured, in general agree- 
ment with expectations. 

It is also interesting to compare to previously known popu- 
lations of Compton thick AGNs. Tozzi et al. (2006) identify 
6 Compton thick AGNs with N H > 10 24 at 1.5 < z < 2.5, de- 
tected by Chandra in the full GOODS-S, which translates into 
a space density at least an order of magnitude lower than that 
inferred here for our mid-IR excess objects. A larger density 
of Compton thick AGNs is directly detected at z ~ 1 by Tozzi 
et al. (2006), but this is still a factor of at least 3 lower than 
our estimate at z = 2 (Fig. fTOb . Overall, the space density 
of Compton thick AGNs that we infer at z = 2 is reasonably 
consistent, within a factor of 2, with the predictions of the 
background synthesis model of Gilli et al. (2007; Fig. ITOb . 
when counting all the Compton thick sources in the models 
with 10 42 < L x < 10 44 erg s -1 , which are found to have an 
average X-ray luminosity of Lx = 1.3 x 10 43 that is similar 
indeed to that of our sources. 4 This large population of Comp- 
ton thick z ~ 2 AGNs has probably similar intrinsic properties 
to that found by Martinez-Sansigre et al. (2006), who select 
z ~ 2 AGNs as mid-IR bright sources with a radio detection 
and a faint IRAC counterpart. They have uncovered a popula- 
tion of AGNs with a sky density of w 10 deg~ 2 , turning out to 
be heavily obscured sources in nearly all cases. Their AGNs 
are likely the tip of the iceberg of the population found here, 
given their sky density some 2 orders of magnitude lower than 
that of our sources, but much higher bolometric luminosities. 
Similarly, Polletta et al. (2006) find a population of luminous 
Compton thick AGNs at z ~ 2 using Spitzer, with a sky den- 
sity of about 25 deg~ 2 . 

In parallel to this work, Fiore et al. (2007) detected a sim- 
ilar population of z ~ 2 obscured AGNs from X-ray stacking 
of galaxies with large 24 pm to optical flux ratios. Given their 
selection criteria, we expect that many of the Fiore et al. ob- 
jects would classify as mid-IR excess galaxies. 

5.3. Contribution to the X-ray background 

Given an observed average hard X-ray flux of ~ 
10~ 16 erg s _1 cm~ 2 (calculated directly from the model 
in Fig. 9 which gives a factor ~ 1.7x higher flux than 
given in Table 1; see Sect. 4.5) and the sky density of 
3200 deg~ 2 , the mid-IR excess sources produce a background 
of 0.7 keV s cm~ 2 sr~' in the 2-8 keV band. This is a quite 
small contribution to the total background in this band, which 
is not surprising given that this band is almost completely re- 
solved by direct detections of individual AGNs (at least up to 
« 6 keV; Worsley et al. 2005; 2006). The relative contribu- 
tion of these mid-IR excess sources increases towards higher 
energies, as these sources are strongly obscured. 

As an example, we have tried to extrapolate the background 
flux emitted in the 10-30 keV band. Using the best fitting 
Gilli et al. (2007) model with N H = 10 24 5 cm" 2 , we esti- 
mate that these sources would produce a background flux of 

4 We emphasize that the space density of Compton thick AGNs at z = 2 
in the Gilli et al. 2007 models is constrained by the shape of the missing 
X-ray background, and is uncertain by a factor of a few. This is because 
z = 2 Compton thick AGN in the models produce only a small fraction of the 
missing background, see Sect. 5.3. 
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about 1.2 keV s _1 cm~ 2 sr~' in the 10-30 keV band. For 
the physically motivated SEDs described in Gilli et al. (2007) 
this prediction turns out to be only weakly dependent on the 
assumed Nn, for AGNs around the Compton thick regime. 
However, for the X-ray spectral shape at the highest energies 
of NGC 1068 and NGC6240 (Vignati et al. 1999), we would 
derive a contribution a factor of 2 larger than above, being 
flatter at the highest energies than the Gilli et al. (2007) mod- 
els (on the other hand, e.g., the Circinus galaxy has a spectral 
shape even softer than the Gilli et al models, Matt et al. 1999). 
These calculations are insensitive to the particular fraction of 
mid-IR excess galaxies which are actually X-ray emitters, as 
the dependence is only upon the stacked X-ray fluxes. 

The total background in the 10-30 keV band is estimated to 
be about 44 keV s _1 cm~ 2 sr _1 . On the basis of the Gilli et al. 
(2007) model, w 70-80% of the 10-30 keV background can 
be ascribed to the sources detected in the deepest X-ray sur- 
veys in the w 0.5-10 keV band. Therefore we conclude that 
this population of z = 2 Compton thick AGNs can account 
for ~ 10 — 25% of the missing background in the 10-30 keV 
band, i.e. of the fraction due to Compton thick AGNs which 
are not individually detected in X-rays. This is in reasonably 
good agreement with Gilli et al. (2007) model calculations for 
the background in this energy range emitted by Compton thick 
AGNs at 1 .4 < z < 2.5. The same models predict that most of 
the rest of the X-ray background is produced by similar AGNs 
at z < 1.4 (although, Steffen et al. 2007 find marginal back- 
ground contributions from X-ray undetected MIPS galaxies 
with 24 /im flux density > 80/iJy at z ~ 1). Our results sug- 
gest that the AGNs in these z < 1 .4 galaxies could be identi- 
fied from the presence of excess mid-IR emission. 

5.4. On the coeval growth of Black Holes and bulges 

About 20-25% of the z ~ 2 star forming galaxies in 
GOODS with K < 22 (i.e. with stellar masses > 10 10 M S ) 
have 24 /im excesses by more than 0.5 dex, with the fraction 
rising to 50-60% for the most massive star forming galax- 
ies with M > 4 x 10 10 M Q (Fig. |5). We have shown that this 
appears to be due to the presence of heavily obscured AGNs 
in these objects, indicating that heavily obscured AGN activ- 
ity was widespread among star forming galaxies and massive 
galaxies in general at z — 2. This is different from the local 
universe, where the fraction of massive galaxies containing an 
AGN is only a few percent (Kauffmann et al. 2004), support- 
ing the idea that the BH growth from the most massive BHs 
must have been higher in the past (see Fig 2 in Heckman et al. 
2004). The AGN activity is likely connected, in some physi- 
cal way, to the intense star formation activity. In fact, locally 
a comparably large fraction of LIRGs and ULIRGs contain 
AGNs with detectable mid-IR excess (Genzel et al 1998; Yun 
et al. 1999 Tran et al 2000; Laurent et al 2000). Therefore, the 
relative fraction of AGNs in LIRGs and ULIRGs might not be 
strongly evolving, while the fraction of AGNs in massive star 
forming galaxies evolves by a substantial factor. Quite obvi- 
ously, this suggests that the epoch of major build-up of stellar 
mass in galaxies at z ~ 2 (e.g., Daddi et al. 2005b; Papovich 
et al. 2006; Franceschini et al. 2006) coincided with a ma- 
jor build-up phase of the central supermassive BHs, as appear 
to be the case in the submillimeter-emitting galaxy popula- 
tion (e.g., Alexander et al. 2005a, Borys et al. 2005). This 
supports the idea that the local correlations between BH and 
galaxy mass were naturally put in place by a concomitant for- 
mation. Is there evidence of such a parallel growth of stellar 
and BH mass in our sample of z ~ 2 galaxies ? 



The rate of BH growth in these objects can be estimated 
from their bolometric luminosities, as Lbol = vMbrc 2 , with 
a typical energy conversion efficiency -q w 0.1 (e.g., Mar- 
coni et al. 2004). Given the large uncertainty in the bolo- 
metric luminosities of the AGNs in mid-IR excess galaxies 
(10 11 ^ Lbol ~ 10 12 Lq) , we infer that the average BH accre- 
tion rate is in the range M B h ~ 0.05-0.57^ J M yr~' (here 
?/o.i = ??/0.1). The mid-IR excess sources are only ~ 1/4 of 
the sample of massive star forming galaxies (i.e., the growth 
in stellar mass is happening in a ~ 4 times larger sample of 
galaxies that remain undetected in hard X-rays). The average 
BH accretion per galaxy may thus be close to 1/4 of the above 
estimate. We are neglecting here the possible additional BH 
growth inside galaxies classified in our simplification as mid- 
IR normal, and BH growth occurring inside individually X- 
ray-detected z — 2 AGNs. The latter term, from individually 
hard X-ray detected AGNs inside z = 2 massive star forming 
galaxies, would increase the BH growth budget by at most a 
factor of 2, and likely much less (see Sect. 4.5). The massive 
star forming galaxies in our full GOODS-S sample have aver- 
age SFR~ 70 M yr _1 (from the UV and soft X-ray). Assum- 
ing explicitly here the more physical IMF of Kroupa (2001) or 
Chabrier (2003), and considering (for comparison to the local 
Universe) that only some 70% of these stars will survive at 
z = 0, we get a ratio of M m /SFR ~ (0.35-3.5) x 10~ 3 %J. 

The local BH to galaxy mass correlation implies 
MBn/M ha \ ge w 1.5 X 10~ 3 (McLure & Dunlop 2001; Marconi 
et al. 2004; Ferrarese et al 2006), within about a factor of 
2 uncertainty. Although uncertain, our estimate of the rela- 
tive growth rates brackets the observed BH to stellar mass ra- 
tio, which tantalizingly suggests that this population of mid- 
IR excess galaxies may indeed represent the major building 
phase of BHs in massive galaxies. It is in fact largely expected 
that massive and star forming z = 2 galaxies are the progen- 
itors of local spheroids and bulges (e.g., Daddi et al 2004ab; 
Shapley et al 2005; Adelberger et al 2005; Kong et al. 2006). 

From the space density of these sources and their BH ac- 
cretion rate, the cosmic BH accretion rate density turns out to 
be 10~ 5 - 10 Al.yr '.\lpc ; . This is a factor of - 2.5-25 
larger than the contribution estimated for z ~ 2 bright SMGs 
(Alexander et al. 2005a), which appear to represent the most 
extreme star-forming galaxies at this epoch; clearly, fainter 
submm observations may reveal significantly more AGNs and 
BH growth and eventually overlap with the Zf-selected pop- 
ulation investigated here (which have S^o^m ~ 1-1.5 mJy, 
see Daddi et al 2005b and Paper I). However, this accretion- 
rate density is comparable to the contribution of luminous 
quasars at these epochs (e.g., Yu & Tremaine 2002; Alexan- 
der et al. 2005a; Croton et al 2006). Integrating over the 
2 Gyr time span corresponding to the 1.4 < z < 2.5 red- 
shift range, this would produce a local BH mass density 
of (0.2-2) x 1O 5 M Mpc~ 3 . Comparing to the estimate of 
4.6 x 1O 5 M Mpc~ 3 for the local integrated BH mass den- 
sity (Marconi et al. 2004), our estimate is about 5-50% of 
the total. This is similar to the estimate that today's massive 
galaxies are forming most of their stars in the 1.4 < z < 2.5 
redshift range, accounting for 1=3 20% of all the stars in the 
local Universe (Daddi et al 2005b). 

5.5. AGN duty cycle 

The high incidence of AGNs inside massive star forming 
galaxies, 20-30% for the overall samples, directly implies that 
the AGN activity at the luminosity levels at which we detect it 
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has a relatively long duty cycle, at most only 3-4 times shorter 
that that of star formation. If only a fraction of mid-IR excess 
galaxies are actually X-ray luminous these estimates would be 
reduced by the same fraction. Accounting for individual X- 
ray detections of massive star forming galaxies would instead 
marginally increase these estimates (Sect. 4.5). 

In Paper I we show that massive star forming galaxies 
with M > 1O U M are generally ultra luminous IR galaxies 
(ULIRGs) with a duty cycle of at least 40% and an estimate 
time duration of at least 400 Myr. At these high masses the 
fraction of mid-IR excess sources reaches ~ 50%. Therefore, 
our results suggest that this obscured AGN activity in massive 
(M > 10 1 1 Mq) galaxies has a duty cycle of at least <~ 20% and 
typical durations of at least 200 Myr. We emphasize that these 
figures are likely to be lower limits because, given the down- 
sizing trends, actively star forming galaxies are likely to rep- 
resent a higher fraction of all massive galaxies at masses lower 
than 10 Mq at z = 2 (i.e., the fraction of quiescent/passive 
galaxies likely increases rapidly with stellar mass at z = 2). 
This is consistent with the timescales inferred by Marconi et 
al (2004). 

It is also worth noting that the strong increase of the frac- 
tion of mid-IR excess galaxies with stellar mass (see Fig. 5) 
might imply that the duty cycle of AGN activity inside star 
forming galaxies depends on stellar mass. This could actually 
be the ultimate reason for downsizing, i.e. the fact that mas- 
sive elliptical galaxies are the first to stop their star formation 
activity (e.g., Cimatti et al. 2006; Bundy et al. 2006; Scarlata 
et al. 2007). With more massive systems being more likely to 
host an AGN, they indeed would be the first ones to turn pas- 
sive if AGN activity is the preferred feedback channel for the 
quenching star formation. Kriek et al (2007) have proposed a 
similar argument. 

5.6. Submm galaxies, merging and BH growth 

It is relevant at this point to summarize the comparison be- 
tween BzK galaxies, particularly the mid-IR excess objects, 
and the co-evally selected SMGs, as this has possibly interest- 
ing implications for the triggering of BH and galaxy growth. 
We have discussed in Sect. 5.4 that the BH accretion-rate den- 
sity at z = 2 is higher by roughly an order of magnitude in the 
mid-IR excess galaxies, over what found in SMGs. Estimates 
of the z = 2 star formation rate density also show that SMGs 
account for a similar share of the overall stellar mass growth. 
Compared to the mid-IR excess galaxies selected to K < 22, 
SMGs are much rarer (at least at bright submm flux densi- 
ties of /sso^m > 5 mJy; e.g., Coppin et al. 2006), have stellar 
masses larger by roughly a factor of 2-3 (e.g., Greve et al 
2005), have star formation rates about 1 order of magnitude 
higher (« 700 M Q yr _1 versus w 70 M© yr _1 ; Chapman et 
al. 2005; this paper) but AGN luminosities are quite similar 
(a few 10 43 erg s _1 ; Alexander et al 2005ab; this paper). This 
can explain why generally the SMGs do not display a mid-IR 
excess over their bolometric luminosity inferred from radio 
(Pope et al. 2006; see also Menendez-Delmestre et al. 2007; 
Valiante et al. 2007), as the mid-IR emission from the AGN 
is likely overshadowed by that from star formation. This also 
implies that the ratio of BH growth to stellar mass growth is 
higher in mid-IR excess galaxies than in SMGs. For the for- 
mer, the concurrent BH and stellar mass growth can account 
for the buildup of galaxies lying on the locally defined black- 
hole-spheroid mass relation, while this is likely not the case 
for the activity of SMGs (Alexander et al. 2005 ab). 

Submm bright galaxies are likely a special phase during the 



buildup of massive galaxies corresponding to both the highest 
SFRs and specific SFRs. They lie outside the stellar mass ver- 
sus SFR correlation, as described in Paper I, in the same way 
that ULIRGs are outliers of the local stellar mass to SFR cor- 
relation. This is consistent with the idea that SMGs like local 
ULIRGs are triggered by major mergers, as inferred also on 
the basis of their morphologies (Conselice et al 2003; Chap- 
man et al 2003), and as suggested also by the compact sizes of 
the starbursts observed in CO lines and in the radio (Tacconi et 
al. 2006; Bouche et al 2007). Instead, the fact that BzK galax- 
ies in general define a tight stellar-mass to SFR relation argues 
against an important role of major mergers in triggering their 
activity. The long duty cycle of their SFR activity (Daddi et 
al 2005b; Paper I) also supports this idea, while SMGs are 
thought to be comparably much shorted lived (e.g., Greve et 
al 2005 but see also Swinbank et al. 2006). 

Therefore, an important conclusion that one might derive 
from our results is that BH growth in galaxies appear more 
directly connected with relatively stable and long lived star 
formation, despite still at considerably high rates, rather than 
being predominantly associated to rare and short lived ma- 
jor mergers events. We notice that most models of coeval 
BH and galaxy growth and of AGN feedback postulate in- 
stead that major mergers are the most relevant events for BH 
growths and for triggering feedback (e.g., di Matteo et al. 
2005; Schawinski et al. 2006; Croton et al. 2006). Our re- 
sult is, in fact, somewhat surprising given that mergers are 
known to be very efficient in channeling the gas toward the 
center of the galaxies. Major mergers are likely still required 
to justify the rapid morphological transition between massive 
high-z star forming galaxies with irregular/disk morphologies 
into bulge dominated systems (although perhaps disk insta- 
bilities might also play this game; Genzel et al 2006), but it 
appears that a substantial fraction of the BH masses inside 
local massive galaxies might have been put in place outside 
such events. 

6. AGN FEEDBACK VIA COMPTON HEATING ? 

The widespread presence of obscured AGNs lurking inside 
massive star-forming galaxies naturally leads to the question 
of their possible feedback effect on their host galaxies. Feed- 
back from AGNs is widely believed to be the key factor for 
solving a major problem of galaxy formation models based 
on ACDM cosmology. In short, these models require an ad 
hoc mechanism to quench star formation in massive galax- 
ies, which is empirically required to reproduce the observed 
color bi-modality of galaxies as well as the old ages of stars 
in early-type galaxies. Passively evolving elliptical galaxies 
are now known to be present in significant numbers at least 
up to z ~ 2.5 (Cimatti et al. 2004; McCarthy et al. 2004; 
Daddi et al. 2005a). As first shown by Granato et al. (2001; 
2004), physically motivated semi-empirical models of galaxy 
evolution implementing AGN feedback can successfully ac- 
count for the existence of massive and passive galaxies al- 
ready at early epochs (see also de Lucia et al. 2006; Springel 
et al. 2006; Croton et al. 2006; Bower et al. 2006; Kitzbich- 
ler & White 2006), as shown also in Paper I. Our results 
are consistent with the idea that potentially all z = 2 mas- 
sive star forming galaxies contain AGNs, and therefore the ad 
hoc mechanism advocated by the models appears at least to 
be plausible because the AGNs are indeed there. Therefore, it 
would be of great interest to find a direct observational con- 
nection between AGN activity and star-formation quenching 
in massive galaxies. 
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We can imagine two ways that energy can be transferred 
from the AGN to the host interstellar medium (ISM), heat- 
ing it up until it escapes from the galaxy in a wind. One is 
via (sub)relativistic jets, and their possible role in a similar 
context - suppression of galaxy cooling flows - has been in- 
voked by, e.g., Binney & Tabor (1995). Alternatively (or in 
addition), the hard X-rays emitted by the AGN may produce 
Compton heating as they are absorbed by the ISM, again pow- 
ering galactic winds (Ciotti & Ostriker 1997; 2007). AGN jets 
and accretion related outflows are likely a major channel for 
mechanical feedback in distant galaxies. However, we cannot 
state anything new on the issue based on our sample of z ~ 2 
galaxies and our results. We focus instead here on the lat- 
ter option, as the likely Compton thick nature of these AGNs 
indicates that Compton heating could indeed be taking place. 

The large obscuration of the central BH implies that very 
high energy photons are heavily absorbed by the material sur- 
rounding the AGN, which can thus be efficiently heated. The 
typical star forming galaxy in our GOODS-S sample has a 
stellar mass of M <~ 3 x 1O 1O M , and the binding energy of the 
gas can be estimated to be of order of M gas <7^ «3x 10 5& f gas , 
where f gas is the gas fraction of the galaxy's baryonic mass 
and <r v is the velocity dispersion. The X-ray luminosity ab- 
sorbed at high energies, e.g., > 2 keV in the rest-frame, is 
10 43 erg s l . These very high energy photons can effi- 
ciently heat the gas, through Compton scattering. Assuming, 
as a limiting case, that all the absorbed energy at > 2 keV goes 
into extracting gas from the galactic potential well, and using 
fgas = 0.1, it would take only ~3x 10 6 yr to blow-off the gas 
completely. Clearly, the efficiency is likely to be much lower, 
chiefly because the material in the BH surroundings can ab- 
sorb and re-radiate to lower energies most of the energy of 
the absorbed hard X-ray photons. The mid-IR excess itself 
testifies for (part of) the AGN energy output being indeed de- 
graded. The energy available for feedback is the difference 
between the total hard X-ray energetic output of the AGN, 
and the fraction of it which is re-radiated at lower energies. 
Since both terms are affected by large uncertainties, estimat- 
ing their difference would be affected by even larger uncer- 
tainties. Notwithstanding this limitation, the large unobscured 
luminosity in hard X-rays that we derive for these objects ar- 
gues for the plausibility of a scenario in which Compton heat- 
ing plays an important role in quenching star formation. 

7. FUTURE PROSPECTS 

We have discovered a substantial population of heavily ob- 
scured, Compton thick AGNs, residing in a large fraction of 
massive star forming galaxies at z ~ 2. We have shown that 
this result has important implications for BH-galaxy evolu- 
tion but our analysis still leaves open many questions. There 
are two crucial points that will need to be addressed by future 
studies: a cleaner selection of mid-IR excess objects and an 
improved estimate of the obscuring column densities, hence 
ofZ4 GN ,andofZ£gN 

For the selection, the main current limitation is that the 
fiducial SFR activity of the host galaxy, used to diagnose the 
presence of an excess at mid-IR bands, is derived from the 
UV. Due to the necessity of large (hence uncertain) redden- 
ing corrections, and in some case of UV obscuration exceed- 
ing what can be inferred from the UV slope, this procedure 
is far from being ideal. A better procedure would be to de- 
rive a reliable measurement of the bolometric IR luminosity 
of the galaxy from long wavelength observations, and com- 
pare this to the mid-IR rest-frame luminosity. The availability 



of ALMA, SCUBA2 and especially Herschel observations in 
the near future should all substantially help with this task, pro- 
viding much better multiwavelength SEDs for the dust emis- 
sion, from mid-IR wavelengths through the far-IR peak and 
beyond. With such a detailed SED, one could try to disentan- 
gle the cold-dust and the warm dust emission, due to, respec- 
tively, star formation and AGN heating. 

Measuring reliable values for Nh, determining accurately 
what fraction of mid-IR excess sources are Compton thick, 
and perhaps also detecting individually these sources, would 
require extremely sensitive observations. With Chandra this 
could require an ultra-deep integration, of order of 5-10 Ms, 
to allow to detect some of the brightest objects. Perspectives 
also look promising for the next generation of X-ray satel- 
lites currently under discussion, such as XEUS (Bavdaz et al 
2006); if realized with a fairly high angular resolution (~ 2- 
3") XEUS would be ideally suited to X-ray spectroscopy of 
these X-ray faint AGNs. Rest-frame optical spectroscopy can 
also provide insights by looking for luminous high-excitation 
emission lines such as [OIII]A5007. 

Mid-IR spectroscopy will also be crucial to improve our 
understanding of these sources. We have not as yet explicitly 
demonstrated that AGN emission does account for the entirety 
of the mid-IR excess flux density. Our analysis cannot exclude 
the possibility that intrinsically different mid-IR SED proper- 
ties, e.g. anomalously strong PAHs, exist in a fraction of z = 2 
galaxies, with the mid-IR still powered by SFR. This might be 
related to the presence of different SFR modes, rare or non- 
existent by z < 1, a result that would be quite interesting on 
its own. Clearly, the ideal means to obtain a direct diagnostic 
of the cause of the 24 ^m excess would be to perform spec- 
troscopy at these wavelengths with Spitzer and Akari. Given 
that the typical mid-IR excess galaxy is highly star forming, 
we do expect to see PAH features in a large majority of the 
mid-IR excess galaxies. However, if the mid-IR excess is di- 
rectly produced by the AGN, we should detect also an impor- 
tant contribution from AGN continuum in most of the mid-IR 
spectra, lowering the equivalent width of the PAH emission 
features. As mid-IR excess objects even at z ~ 2 are relatively 
bright at 24 ^m, with flux densities often in excess of 100- 
200^Jy, the observations are feasible and useful spectra are 
already available in the Spitzer archive, although their analy- 
sis exceeds the scope of the present paper. 

8. SUMMARY AND CONCLUSIONS 

We have studied the nature of a population of z ~ 2 galax- 
ies that displays a marked excess flux density at mid-IR wave- 
lengths, over what expected to arise from just normal star for- 
mation activity as estimated from their rest-frame UV. The 
examined sample includes over 600 /f-selected z ~ 2 galaxies 
in the two GOODS fields (with 25% of them having a spectro- 
scopic redshifts), with multiwavelength informations includ- 
ing data in UV, optical, near-IR, Spitzer IRAC (3.6-8.0/im 
bands), Spitzer MIPS (24 fim and 70 ^m bands), SCUBA, 
VLA 1.4 Ghz, and X-ray from Chandra. This sample has 
been cleaned by all known AGNs, identified via their Chandra 
2-8 keV X-ray emission and/or by the presence of a power- 
law SED over the IRAC bands. The major results of this work 
can be summarized as follow. 

• Roughly 20-30% of the massive star forming galaxies 
at z ~ 2, selected down to K < 22, display a mid-IR ex- 
cess larger than about a factor of 3, which we classify 
as "mid-IR excess galaxies". This fraction increases 
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rapidly with the stellar mass of the galaxies, reaching 
- 50-60% for M > 4 x 1O 1O M . We have shown that 
our results do not strongly depend on the exact thresh- 
old used for defining these mid-IR excess galaxies. 

This population of mid-IR excess galaxies is at slightly 
higher average redshift than normal galaxies (<z>— 
2.1 vs. 1.9), but overall the two redshift distribu- 
tions are similar. This population becomes domi- 
nant at the highest 8 /jm (rest-frame) luminosities (for 
L(8^m) > 10 L Q ), that we accurately derive at z ~ 2 
from Spitzer MIPS 24 /mi flux densities. Relatively 
shallow survey with Spitzer MIPS of the z ~ 2 Universe 
are thus the ones that are proportionally most affected 
by the presence of mid-IR excess galaxies. 

Mid-IR excess sources have redder K — 5.8 colors than 
normal galaxies, reflecting different intrinsic SEDs in 
the optical to near-IR rest-frame, that can not be at- 
tributed to mere redshift effects. For these mid-IR ex- 
cess sources, the peak of the stellar dominated part 
of the SED is longward of the canonical 1.6/tm rest 
frame. We interpret this as evidence for the presence 
of a warm dust AGN contribution to the continuum in 
these sources, consistent also with the finding that the 
excess flux density seems to become weaker beyond 
20/im rest-frame, based on Spitzer 70 /im observations. 

This excess AGN component over the near-IR rest- 
frame will likely bias the estimates of stellar masses and 
photometric redshifts for some of the z ~ 2 galaxies, 
when including IRAC photometry in the SEDs. How- 
ever, we quantify the effect on the derived photometric 
redshifts to be only at the level of 2-3% in (1 +z), for 
the mid-IR excess sources. 

Stacking of Chandra X-ray data reveals a strong di- 
chotomy in the two populations: the normal sources 
have a soft spectrum, consistent with star formation in 
both 0.5-2 and 2-8 keV bands. Contrary to them, mid- 
IR excess objects exhibit a much harder spectrum, un- 
ambiguously showing the presence of highly obscured 
AGNs inside these sources. 

We use a variety of methods to infer Nh, the col- 
umn density of the absorbing material along the line 
of sight. Our best estimate suggests typical obscuration 
of A^h = 10 24 - 5 cm~ 2 , derived from the X-ray spectrum 
after accounting for the X-ray emission due to star for- 
mation. As a result, we infer that the majority of the 
mid-IR excess galaxies is likely to host a Compton thick 
AGN. 

This previously unknown population of AGNs has typ- 
ical X-ray luminosity of (1-4) x 10 43 erg s _1 in the 2- 
8 keV rest frame band. We infer a range for the bolo- 
metric luminosity of 10 11 — 10 12 L Q . This is comparable 
to the bolometric luminosities of their host galaxies, as 
due to star formation, for which we estimate a typical 
value of -5 x 10 n L©. 

The sky density of Compton thick AGNs inside massive 
star forming galaxies to K < 22 is ~ 3200 deg~ 2 and 
their space density is <~ 2.6 x 10~ 4 Mpc~ 3 . Their space 
density is higher than that of all AGNs already known 
at z = 2, and much higher than that of previously known 



Compton thick populations at high redshifts. The space 
density of Compton thick AGNs that we derive agrees 
reasonably well with that predicted by the background 
synthesis models of Gilli et al (2007). 

• Despite their large space density, we extrapolate that 
this population can account for only about 10-25% 
of the missing X-ray background in the observed 10- 
30 keV band, in agreement again with the prediction of 
Gilli et al (2007). Most of the high energy background 
still to be resolved is likely inside similar sources but 
residing at z ^ 1.4 rather than at z ~ 2. 

• The widespread presence of obscured AGNs inside 
massive star forming galaxies at z ~ 2 seems to quite 
naturally account for the concurrent growth of the BH 
together with the stellar mass of the hosting galaxies. 
Within quite large quantitative uncertainties, this pic- 
ture is consistent with these objects approaching the lo- 
cal values for the ratios between the mass of central BH 
and that of their hosting bulges. 

• A relatively long duty cycle for the AGN activity at 
z = 2, only a factor of 3-4 lower than that of star forma- 
tion activity, is suggested by the widespread presence 
of AGNs that appear to be present among massive star 
forming galaxies. 

• Comparing to submm selected galaxies (SMGs), we 
find that BHs are growing faster in mid-IR excess 
galaxies, relatively to the ongoing SFR, and their 
integrated contribution to the BH accretion-rate density 
are also larger. This suggest that of BH growth in 
massive galaxies is taking place outside major merger 
events. 



• The large, increasing incidence of Compton thick 
AGNs with stellar mass, along with the energy deposi- 
tion provided by the absorption of the hard X-ray pho- 
tons, suggest that these objects may be responsible for 
the AGN feedback eventually leading to the termina- 
tion of the star formation activity, and to do so begin- 
ning with the most massive galaxies as in the downsiz- 
ing scenario. Although existing data do not allow to 
firmly prove this tantalizing possibility, we show that 
the involved energetics is consistent with this picture. 

• We finally discuss future prospects for improving our 
still limited understanding of the properties of mid-IR 
excess galaxies and Compton thick AGNs at high red- 
shift. The crucial and obvious required next step is the 
full construction of their spectral energy distribution, 
including ALMA, SCUBA-2, and Herschel facilities as 
they become available, as well as of their spectral mid- 
IR properties, via Spitzer IRS and Akari spectroscopy. 
Extremely deep Chandra (5-10 Ms) exposures would 
be required to directly detect some of these Compton 
thick AGNs. 

In summary, we have uncovered a major population of lu- 
minous and obscured z = 2 AGNs, previously largely un- 
known. With AGN and star formation activity in massive 
galaxies peaking at z ~ 2, this newly detected population of 
Compton thick AGNs appears to be an excellent candidate 
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for playing a major role both for establishing the BH mass- 
bulge mass relation, and for providing the feedback necessary 
to discontinue further star formation in massive galaxies. 
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